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Abstract 

A method based on a differential algebraic equation (DAE) approach is em¬ 
ployed to find stationary solutions of the Poisson-Nernst-Planck equations with 
steric effects (PNP-steric equations). Under appropriates boundary conditions, the 
equivalence of the PNP-steric equations and a corresponding system of DAEs is 
shown. Solving this system of DAEs leads to the existence of stationary solutions 
of PNP-steric equations. We show that for suitable range of the parameters, the 
steric effect can produces infinitely many discontinuous stationary solutions. More¬ 
over under a stronger intra-species steric effect, we prove that a smooth solution 
converges to a constant stationary solution. 


1 Introduction and statement of results 


To study biological ion channels, the Poisson-Nernst-Planck system (0, a, 0,0) is 
commonly used as a model to describe the ionic flows in open ion channels. The Poisson- 
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Nernst-Planck system is given by 
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for X e (— 1 , 1 ) and t > 0 . Here Cn = Cn{x,t) is the density of anions and Cp = Cp{x,t) 
is the density of cations; 0 = 0(x, t) is the electrostatic potential; Dn and Dp are the 
diffusion coefficients; —Zn and Zp are positive integers; p is the permanent charge density; 
£ is a parameter. Next we introduce a new type of Poisson-Nernst-Planck system by 
considering the steric effect (or size effect) which occurs due to the fact that each atom 
within a molecule occupies a certain amount of space [lailDlElIISlIls]: 
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( 1 . 2 ) 

for X G (— 1 , 1 ) and f > 0 . Here c„ = Cn{x,t) is the density of anions and Cp = Cp{x,t) is 
the density of cations; 0 = (j){x, t) is the electrostatic potential; A = A{x) is the cross- 
sectional area of the ion channel at position x; Qnn, Qpp and Pnp are positive constants; 
—Zn and Zp are positive integers; 5 and e are two parameters. 

Assuming constant cross-sectional area A(x), without loss of generality, we consider 
the following two-component drift-dijfusion system 


' Ut = V ■ {di Vu + 'diu V 0 ) -|- V ■ [pii u Vu -I- gi2 u Vv), x G H, t > 0 , 

< Vf = V ■ (d2Vv -h 'd2'cV(/>) + V ■ (g2ivVu -h g22'cVv), x G hi, t > 0 , ( 1 . 3 ) 

A 0 = 7iM-I-72X, X G n, t > 0 , 

where u = u{x,t) and v = v{x,t) are densities of the two species u and v, which are 
assumed to be nonnegative functions; 0 = 0(x, t) is the electric potential; di and ^2 are 
diffusion rates. Throughout this paper, "di, ^fn, gi2, g2i, g22 and 71 are positive constants; 
'§2 and 72 are negative constants, unless otherwise specified, hi C M"' is a bounded domain 
with smooth boundary. 

When the domain under consideration is extended to the entire space M”, fll.dp be¬ 
comes 

' Ut = V ■ {di Vu A diu V 0 ) -|- V ■ {gn u Vu + gi2 u Vx), x G M", t > 0 , 

< Xt = V ■ (^2 Vx 52 X V 0 )V ■ (5'2i "y Vm 5^22 "y Vx), X G M”, t > 0 , ( 1 . 4 ) 

A 0 = 7i X-I-72X, X G M"', t > 0 . 
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We note that, in the absence of guuVu, gi^uVv, g2iuVv, and (722 'fVn, fll.3p reduces 
to the Keller-Segel type equations 

"wt = V • (di Vn +'di M V(/)), x e M", t > 0 , 

< Ut = V ■ (^2 Vn + ^2 n V(/)), a; e M"', t > 0 , ( 1 . 5 ) 

^ — A 0 = 7 i M + 72 n, X e M"", t > 0 . 

Keller-Segel model is a classical model in chemotaxis introduced by Keller and Segel 
([E])- For the last two decades, there has been considerable literature devoted to the 
Keller-Segel model. For the Keller-Segel model, we refer to the textbook [moK] , review 
papers mm and references therein. 

In this paper, we are concerned with stationary solutions to fll.3p and fll.4p . i.e. with 
time-independent solutions to the following elliptic systems 

"0 = V ■ (di Vu + 'diu V 0 ) -I- V ■ {gii u Vu + gi2u'Vv), x e Q, 

< 0 = V ■ ((^2 Vn-I-^2 n V 0 )-I-V • ((721 "y Vu-I-(722 "y Vn), x e fl, ( 1 . 6 ) 

^ — A0 = 7 i M 4- 72 n, X e fl. 


Using the elementary fact that V(logM) = Vn/n, the hrst and second equations in 
m can be rewritten as 


{ 0 = V • (uV{di \ogu + -di (j) + giiu + gi 2 v)j, x e fl, 

0 = V ■ V(d 2 logn -4 ^2 0 + 5^21 M + 5'22 "y) j, x e fl. 

It is readily seen that if we can hnd u, v and 0 satisfying the algebraic equations 


(1.7) 


di\ogu + 'di(f) +giiu +gi2V = Cl, x e fl, 
d2 logn 4 - ^2 0 + g2i u + g22V = C2, x e fl. 


( 1 . 8 ) 


where Ci and C 2 are constants, then such u, v and 0 automatically form a solution of 
fll.7p . A natural question arises as to whether any solution of fll.7p also satishes fll.8p . 
It will be shown in Proposition 12.11 that the answer is indeed affirmative when certain 
appropriate boundary conditions are imposed on the solutions, i.e. 


and 


where Fi = 
instance, ^ 


u Fi 


dFi 

dv 


< 0 , 


on dVL, 


(1.9) 


nF2-;:^<0, on dVL, (1-10) 

ou 

di \ogu + 'di4> + gii u + gi 2 V and F 2 = ^2 logn -4 ?92 0 + 5'2i u -4 g 22 v. For 
= 0 (i = 1, 2) on dQ, or the Neumann boundary conditions 


du dv 
do du 


(90 

du 


0 on dQ. 


( 1 . 11 ) 
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As a consequence, our problem now turns to establishing the existence of solutions 
for the differential algebraic equations (DAEs): 

'di\ogu + i!}i(p + gnu + 9 i 2V = Cl, x e ff, 

< d2\ogv + 'd2 4> + 92iU +g22V = C2, X e n, (1-12) 


= 7iM + 72^, xeQ. 

To the best of the authors’ knowledge, this is the hrst time such DAE approach has 
been employed to drift-diffusion systems such as fll.hp . In Section [21 we investigate in 
Theorem 12.21 the existence and uniqueness of solutions u = m( 0) and v = n(0) to fll.Sp 
under the following hypothesis: 

(HI) 911 922 — 912 921 > 0. 

Note that fll.Sp is a system of nonlinear algebraic equations for which explicit solutions 
expressed by the form u = u{(j)) and v = v{(j)) in general cannot exist. Due to (HI) 
however, the solution u = u{(j)) and v = n(0) of fll.81) in implicit form can be given 
uniquely. With the aid of Theorem 12.21 we are hnally led to the following semilinear 
Poisson equation 

—A0 = G{(p), X G D, (1-13) 

where G(0) = 7 im( 0) -|- 72n(0). To establish existence of solutions of fll.lhp under the 
zero Neumann boundary condition 

^ = 0 on 5D, (1-14) 

ou 

more delicate properties of the nonlinearity G = G{(j)) and the solution u = «(</>), v = v{(j)) 
are explored in Proposition 12.31 These important properties turn out to be essential in 
proving the existence of solutions of fll.lhp by the standard direct method in the calculus 
of variations. 


Theorem 1.1 (Existence of stationary solutions to PNP-steric equations under 
(HI)). Assume that (HI) holds and ci and C 2 be fixed. (USD coupled with the Neumann 
boundary conditions (11.lip has a solution {u{x),v{x),(j){x)) G C'^(D). 


The intuition behind the differential algebraic equation approach we use in obtaining 
Theorem 1 1.1 1 is elementary. However, the result is remarkable in that only (HI) is needed 
to ensure the existence of solutions to the elliptic system (II. 6p together with the Neumann 
boundary conditions (II.lip . On the other hand, under the hypothesis: 

(H2) 911 922 — 912 921 < 0, 

(II. 8p admits triple solutions (see Section [2D- Moreover, in this case there may exist 
inhnitely many solutions for (II.6p . 
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Theorem 1.2 (Existence of stationary solutions to PNP-steric equations under 

(H2)). Assume that (H2) holds and ci and C 2 be fixed. (USD coupled with the Neumann 
boundary conditions fll.lljl has either a solution or infinitely many discontinuous 
solutions (u(x), f (x), 0(x)). 


2 Two species equations 


To begin with, we show that under the boundary conditions (11.9p and (ll.lOp . every 
solution of (II.7p also solves (II.8p . as mentioned in Introduction. 


Proposition 2.1 (Equivalence of algebraic equations and differential equations). 

The systems of PDEs 

0 = V • {uV {di \ogu + di + Qiiu + gi 2 v'^, x e Q, 


0 = V 


^nV(d 2 logn + 'd 2 0 + 5'2iM + 5'22w)j, x e fl, 


together with the boundary conditions 


and 


uFi ——< 0, on oil, 
ah' 


dFo 

vF 2 ——< 0, on dTl, 


du 


where Fi = di logw + Aifi + guU + gi2 v and F2 = ^2 logn + ■^2 0 + (?2i w + 5^22 v, is 
equivalent to the system of algebraic equations 

di log M + ■di 0 + S'!! M + 5'12 "y = Cl, X G fl, 
d 2 log n + ■^2 0 + 5'21 U + g 22 V = C 2 , X e fl. 

Proof. Integration by parts leads to the desired result. Indeed, 

[ u FiAFi dx = — [ V(m Fi) ■ VFi dx + [ uFi - 9 - ds 
Jn Jn Jan on 

< — / u\VFi\^dx— / FiVu-VFidx. 

Jq Jn 

( 2 . 1 ) 


However, V ■ (m VFi)= 0 gives 0 = Fi V ■ (u VFi) = Fi Vu ■ VFi + u FiAFi. This shows 
that J^u iVFipdx < 0, and thus Fi should be a constant independent of x. In a similar 
manner, we can prove that F 2 is a constant independent of x from V ■ {v VF 2 ) = 0. The 
proof of the converse is trivial. Hence the proof is finished. □ 
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Since we have established Proposition 12.11 it is now necessary to investigate existence 
and uniqneness of solutions to fll.Sjl . In particular, the solution of interest takes the form 
{u,v) = (M( 0 ),n( 0 )). 


Remark 2.1. When the zero Neumann boundary conditions i^ = i^ = i^ = 0on dVL are 
considered, it is easy to see that these boundary conditions lead to ^ = 0 and ^ = 0 
on dVL. However, note that the converse is not true. Since the hrst two equations of fll.3p 
are in divergence form, the total charges are conserved: 


u = 


u{x,t)dx= / u{x, 0 )dx, 


V = 


v{x,t)dx= / v{x, 0 )dx, 


( 2 . 2 ) 


where {u{x,t),v{x,t)) is a solution of fll.3lh It follows from the divergence theorem and 
the Poisson equation in fll.3p that in the case of the homogeneous Neumann boundary 

condition ^ = 0 , we have 

7 ih = - 72 n, (2.3) 

which is equivalent to the electroneutrality condition. Here 71 and 72 stand for the valence 
numbers of ions u and v, respectively. 


Theorem 2.2 (Existence of solutions to fll.Sp b Assume (HI). Then for any given 
0 = 00 G there exists a unique solution {u,v) = (uo,vo) of 


di \ogu + 'di(j) + gnu + guv = ci. 


d2 log n + ■^2 0 + 5'21 u + g22V = C2. 


(2.4) 


Moreover, fl3.8p has a unique solution {u, v, 0) which can he represented implicitly as 
{u,v) = (M(0),n(0)) and M(0),n(0) are functions. Prove that 

Proof. For any given 0o G M, we first show that the system of algebraic equations fl3.8p 
has at least one solution (u, v) = {uq, Vq). In other words, the two curves di \ogu + giiu + 
912 w = Cl — 'di 00 and ^2 log v + g 2 iu + 922 ^ = C 2 — ^2 0o on the v-u plane has at least one 
intersection point. To see this, we differentiate di logM(n) + 911 u{v) + 912 V = ci — Ai 0o 
with respect to v to obtain 


u\v) 


9i2u{v) 
di + 9iiu{v) 


(2.5) 


which is always strictly less than zero because of the term logu appearing in the hrst 
equation of 03.81) and u cannot take nonpositive value in logu. Therefore, the prohle 
u = u{v) of di logu + 9 iiu + 912 V = Cl — di 00 on the v-u plane has the following 
property: 


• as n — 00 , M —)■ 00 ; 

• as n — 00 , M ^ O’*"; 
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• u = u{v) is decreasing in n G M. 


In a similar manner, we have for the second eqnation ^2 logn(M) + g 2 iu + 922 v(u) = 
C 2 - '&2 00 of fl3^ . 


v\u) 


921V 

d2 + 922 V 


< 0 . 


( 2 . 6 ) 


The prohle v = v{u) of ^2 logn + 921 U + 922 v = C 2 — 1^2 (t>o on the v-u plane enjoys the 
following property: 


• as M —— 00 , n —)■ 00 ; 

• as u ^ 00 , V ^ 0+; 

• V = v{u) is decreasing in n G M. 

As a conseqnence, it follows from the property of the prohles of the two curves 
di hgu + 9 nu + 912 v = Ci - ■di 0o and ^2 logn + 921 u + 922 v = C 2 - '^2 0o, that the 
two curves in the hrst quadrant of the v-u plane intersect at least once. That is, given 
any 0 = 0o £ 1^, we can hnd at least one solution {u,v) = {uo,vo) which satishes fl3.8p . 
We eliminate the possibility of non-uniqueness of solutions {u, v) to fl3.8p for a given 
0 G M by contradiction. Suppose that, contrary to our claim, there exist in the hrst 
quadrant of the v-u plane two distinct solutions (Mi,ni) and (M 2 W 2 ) which satisfy fl3.8p 
for given 0 = 0o G M. Denote by Mi{u,v) (respectively, M 2 {u,v)) the slope of the curve 
di log u + 9 uu + 912 V = Cl - di 00 (respectively, ^2 log v 921 U + 922 V = C 2 - '&2 0 o) at 
{u,v). It is easy to observe that. 


(Mi(mi, Ui) - M 2 {ui,vi)) (Mi(m 2 , V 2 ) - M 2 {u 2 , V 2 )) < 0. (2.7) 


Without loss of generality, we may assume that Mi{ui, Vi)—M 2 {ui, Ui) < 0 and Mi{u 2 , V 2 ) — 
M 2 {u 2 ,V 2 ) > 0. According to the Intermediate Value Theorem, there exists {u*,v*) for 
which Mi{u*, V*) — M2{u*, v*) = 0, where u* lies between ui and U 2 while v* lies between 
Vi and V 2 . However, using fl2.5p and (12.61) . Mi{u*,v*) = M 2 {u*,v*) leads to 

912 u* ^ d2 + 922 V* ^2 8) 

di + 9 n u* 921 V* 

It turns out that the last equation is equivalent to 




(2.9) 


which contradicts (HI). Consequently, for given 0 G M, uniqueness of solutions to fl3.8p 
follows. By applying the local implicit function theorem at each point (m( 0), n(0), 0) 
which satishes fl3.8p . we obtain smoothness of the solution {u,v) = (M(0),n(0)). The 
proof is completed. 

□ 


Certain important properties of solutions to (11.8p are investigated in the next propo¬ 
sition. These properties include non-simultaneousness of vanishing u and v, monotonicity 
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oi u = u{(j)) and v = v{(j)), and monotonicity and convexity oi u = u{v). An important 
consequence of these properties is fl2.15p . which turns out to play a crucial role in em¬ 
ploying the direct method to establish existence of solutions for the semilinear Poisson 
equation fll.131) . 


Proposition 2.3 (Properties of solutions to fll.Sp h As in Theorem 12.21 assume 
(HI). Then fl3.8p is uniquely solvable by the implicit functions {u,v) = {u{(j)),v{(j))). 
Concerning the properties of {u,v) = (m( 0), n(</>)), we have 

(i) (Non-simultaneous vanishing of u = u{(j)) and v = v{(j))). For fixed 0 = 0o 
and {u,v) which solves fl3.8p . for 5 > 0 one of the following holds: 


u > 6 , 

V > 6, 

( 2 . 10 ) 

u > S, 

V < 6, 

( 2 , 11 ) 

u < S, 

V > 6. 

( 2 . 12 ) 


It is equivalent to say: there exists 6 > 0 so that u < 6 and v < 6 cannot simulta¬ 
neously be true; 


(ii) (Monotonicity of u = u{4>) and v = v{4>)) u'{4>) and can be expressed in 
terms of u{(f)) andv^f), i.e. 


u'if) 


gi2 O2 + 922) 

912 921 - + 9ll^ + 922^ 


(2.13) 


v'if) 


921 Oi - 62 + 9ii) 

fi'12 921 - + S'!!) + ^'22) 


(2.14) 


Also, u'{(f)) < 0 and n'(0) > 0 for 0 G M, i.e. u = m(0) is monotonically decreasing 
m 0 G M, while v = n(0) is monotonically increasing in 0 G M. Furthermore, u'{4>) 
and n'(0) are uniformly bounded for 0 G M. In addition, we have 


7in'(0) -F72n'(0) < -k, 


(2.15) 


for some constant k > 0 ; 


(Hi) (Monotonicity and convexity of u = u{v)) u can be expressed by u = u{v). 
Moreover, 


,, ^ _ u{v) {d2 Oi + 922 Oiv - 912 O2 v) „ 
V {di 62 +911(^2 u{v) - g2i9iu{v)) ^ 
which implies competition between u and v, and 

di 62 u'{vf — d2 61 u‘^{v) 


(2.16) 


u"{v) = 


u{v) {di 62 + 9 ii 6*2 u{v) - g 2 i Oi u{v)) 


> 0 . 


(2.17) 















Proof. We prove {i) by contradiction. Assume to the contrary that there exists 5 > 0 
such that u < 6 and v < 6 hold simultaneously. As 6 is sufficiently small, n, v are even 
smaller and log n, log n < 0. On the other hand, Oicf and 6*2 0 have opposite signs since 
6*1 > 0 and 6*2 < 0. Therefore, in fl3.8l) either di logu cannot be balanced with 6*1 0 or 
d 2 log V cannot be balanced with 62 0- Either case leads to a contradiction. 

To prove (ii), we hrst differentiate the two equations in fl3.8p one by one with respect 
to 0 , and obtain two equations in which the unknowns can be viewed as m'( 0 ) 

Solving them gives m'( 0) and n'(0) as stated in (ii). Due to (HI), it is easy to see that 
m'( 0 ) < 0 and n'( 0 ) > 0 for 0 e M. 

Using (i), when u,v > 6 , clearly it follows from 02.131) and 02.141) that 


m'(0) < 


gi2 6*2 ~ (^ + 922) 


912 6*2 — 922 6*1 


and 


9i2 921 ~ (t + 9ii) (t + 922) (x + 9ii) (t + 922) — 9i2 921 

921 6*1 — 6*2 (^ + 9ii) _ 921 6*1 — 6*2 (^ + 9u) 


n'( 0 ) > 


(2.18) 


(2.19) 


9i2 921 - (^ + 9ii) + 922) 9ii 922 - 912 921 

for 0 G M. As for the other two cases u > 6 ,v < S and u < 6 , v > S, it suffices to consider 
one of them since they are symmetric. Suppose that u > S,v < 6 . For u > S,v < Sq 6 , 
02.13P and 02.14p lead to 


lim m'( 0 ) = 


0 i 


0i 


{i^)+ 9 ii) if+ 911) 


and 


lim U(0) = 0 . 


( 2 . 20 ) 


( 2 . 21 ) 


On the other hand, when u > S, 6 o < v < 6 , estimates of m'( 0) and n'(0) similar to 02.18P 
and 02.19P are given by 


m'(0) < 


912 6*2 - 6*1 (^ + 922) 


9i2 6*2 — 922 6*1 


912 921 - (f + 9ii) (1 + 922) (f + 9ii) (1 + 922) - 912 921 


and 


U(0) > 


921 9i — 6*2 (^ + 9n) 


921 9i — O 2 {^ + 9ii) 


912 921 - (^ + 9ii) if + 922) 9ii if + 922) - 912 921 


( 2 . 22 ) 


(2.23) 


As a consequence, we conclude that there exists a constant fc > 0 such thtat 7 im'( 0 ) + 
72 ^'( 0 ) < - k . 

Multiplying the hrst equation in 03. 8 p by 62 and the second equation in 03.8p by 6 * 1 , 
we obtain two equations. Subtraction of the two resulting equations gives an equation in 
terms of u and v. Implicit differentiation then gives the desired result in {in). 

□ 
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Theorem 12.21 and Proposition 12.31 lead us to consider the Neumann problem for the 
semilinear Poisson equation fll.l3jl . i.e. 


{ —A0 = G{(p) in n, 
^ = 0 on dQ. 

We are now in a position to prove Theorem 11.11 


(2.24) 


Proof. (Theorem II.ip 

First of all, it can be shown by means of Proposition 12.31 that G{(j)) satishes (Al) and 
(A 2 ) of Theorem 13.31 Upon using Theorems 12.21 and 13.31 we establish the existence of 
fll.l2p . Applying Proposition 12.11 completes the proof of Theorem 11.11 

□ 


3 Triple solutions; ^^j-induced triple solutions 


For the case where the hypothesis (H 2 ) is assumed we can also establish an existence 
theorem for 01.121) . In particular, in this case 01.Sp admits a triple solution in the sense 
that for some (fo, we can hnd three pairs of solutions (Mj(0o),nj(0o)) (* = 1)2,3). See 
Figure im for an example. This is essentially different from the case where (HI) holds. 




1.0 

0.8 

0.6 

0.4 

0.2 



1.50 1.55 1.60 

Figure 3.1: A triple solution. 


1.65 


<t> 
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As shown in Figure IHTTl there may exist a triple solution and thus uniqueness no longer 
holds when (HI) is violated, i.e. when ^fn g22—gi2 5'2i < 0 . When non-uniqueness occurs, 
it is readily seen that there must be a 0 * G M such that u'{(j)*) = v'{(j)*) = oo. Due to 

and fl 2 . 14 D . this can occur only when gu fi'21- + fi'22) = 0 . Indeed, 

it is possible that the numerator gi2 62-61 5(22) ^ -00 {g2i 61-62 gn^ -)• 

—00) in fl 2 . 13 p f fl 2 . 14 p i as n( 0 ) —)■ 0 + («(</>) —t O’*'). However, as we have mentioned in 
the proof of Proposition 12.31 that 


lim u'{(j)) 

v^0+ 


61 


lim v'{(j)) 

u—>- 0 + 


62 

+ ^ 722 ) 


it is therefore necessary to hnd [u, v) at which gi2 g2i — + 9iij + ^722 j vanishes 

by solving the equations: 


'di logu + '6i4> + giiu + guV = D, 

< ^2 logn + 'd 2 0 + 5 ' 2 iM + 5 ' 22 'y = C 2 , ( 3 . 1 ) 

,5'12 5'21 ~ (^ + S'!!) + 922) = 0 . 

It turns out that fl 3 .ll) is reduced to a single equation a{u) = 0 by eliminating 0 and 

substituting v = f'^i+gn D -_ More precisely, 

o u{gi2g2i—giig22)—dig22 ^ ’ 


. . 1 / , , ^2 (di -|- m) \ 

<^(“) =7- Cl - cfi logu - guU- gi2 -r-- 

6\ \ u [gi2 g2i — gii g22) — di g22 J 


62 


d2 log 


u 


d2 {di + gn u) 

{gi2 fi'21 — gii fi'22) — di g22 ) 

d2 {di + gii u) \ 

- 7721 w — 7722 —^ 1 - 

u [gi2 5'21 — S'!! g22) — di g22 J 


Now the question remains to determine the prohle of cr{u) = 0 . To this end, we hrst 

observe that a(u) makes sense only when u > u* , where u* :=- dL92i - > q _ Also, it 

is readily verihed that 

lim a{u) = lim a{u) = —00. ( 3 . 2 ) 

u^{u*)+ u^oo 

To hnd extreme points of cr(M) = 0 , we calculate 


a'{u) 


_ p{u) (m (^'21 61 - gii 62) - di 62) _ 

61 62 u (di -|- u) {u {gi2 7721 — 77 ii 7722) ~ di 7722)^ 


( 3 . 3 ) 
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where p{u) = + k 2 u‘^ + kiu + ko, and 

^3 = dll {dll 922 — 912 92 iY , 


k 2 — di {gi 2 921 — 3 ^fn g 22 ) {912 921 — 9ii 922 ), 


ki — —di {d 2 921 9i2 + 2 di 921922 9i2 — Sdi gn 5 ^ 22 ) > 


ko = 


dlgl2- 


(3.4) 


We remark that the denominator of <j'{u) is always away from 0 since u > u*. Let the 
numerator of fl3.3p be zero. Then the four roots are 0 > U 2 , and M 3 , where 

Ml, M 2 , and M 3 are the three roots of p{u) = 0. Indeed, mi, M 2 , and M 3 are three distinct 

real roots. To show this, we use Fan Shengjin’s method. As in [ 8 ], let 

A = kl — Skiks, B = kik2 — Qkoks, C = k\ — ?> ko k2 (3.5) 


and the discriminant 

Ai„ = B‘^-4AC. 


(3.6) 


Lemma 3.1 (Shengjin’s discriminant( |I5] )) • There are three possible cases using the 
discriminant Adis ■' 

(i) If Adis > 0 , then p{u) = 0 has one real root and two nonreal complex conjugate 
roots. 

(a) If Adis = 0 , then p{u) = 0 has three real roots with one root which is at least of 
multiplicity 2. 

(Hi) If Adis < 0 , then p{u) = 0 has three distinct real roots. 

It can be shown that Adis A 0 cannot be true under the assumption (H2). Using 
Lemma [ 3 .11 we conclude that the cubic equation 

p{u) = kou^ + k 2 U^ + kiu + ko = Q ( 3 . 7 ) 

has three distinct real roots M3 < M2 < mi. Due to p{±oo) = ±cxo and ko,ko > 0, it 
is easy to see that either M3 < M2 < mi < 0 or M3 < 0 < M2 < mi. However, the case 
M3 < M2 < Ml < 0 cannot occur since a{u) makes sense only for m > m* > 0. For the 
other case M3 < 0 < M2 < Mi, M3 cannot be a extreme point of a{u) because of M3 < 0. 
Accordingly, there are at most two extreme points mi and M2. We have by fl 3 . 2 p the 
asymptotic behavior a{u) —?■ —00 as m —?■ m* or m —?■ 00 , which leads to the fact that 
the number of extreme points of cr(M) = 0 can only be odd. As a consequence, one of mi 
and M2 cannot be a extreme point of ct(m) = 0 and the other one is the extreme point of 
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a{u) = 0. Suppose that U2 is a extreme point of a{u) = 0, then so is ui, which yields a 
contradiction. Therefore, ui is the extreme point of (t{u) = 0. In fact, it can be shown 
that U 2 < u* so that a{u 2 ) dose not make sense. Now a question remains, i.e., how to 
determine the sign of cr(ni)? We see that the maximum of cr{u) is attained at u = Ui or 
maXu^u* (j{u) = (t{ui). Moreover, the criterion for determining the roots of the equation 
a{u) = 0 is stated as: 

• when a{ui) > 0 , a{u) = 0 has two distinct positive solutions; 

• when a{ui) < 0 , a{u) = 0 has no solutions; 

• when a{ui) = 0 , a{u) = 0 has a unique positive solution (i.e. u = ui). 

Noting that as cr(Mi) = 0 and cr'(Mi) = 0, we have m'( 0) = —n'(0) = m"( 0) = —n"(0) = oo 
evaluated at u = ui. In other words u{(j)) and v{(j)) has a reflection point at 0 = 0 for 
some 0 G M. 

On the other hand, as < 0, the denominator gi 2 g 2 i — + h 22 ^ 

of fl2.13p (also fl2.14p l always keeps it sign since a(u) = 0 has no solutions. Indeed, it 
is easy to see that gi 2 g 2 i - + S'!!) + ^' 22 ) < 0 , which leads to m'( 0 ) < 0 and 

n'(0) > 0. We arrive at 


Theorem 3.2 (Trichotomy under (H2)). Assume that (H2) holds and (M(0),n(0)) 
solves 


di \ogu + -dicf) + giiu + gi 2 V = Ci, 


( 3 . 8 ) 


d2 logn + -^2 (/> + 5'21 u + g22V = C2. 


Then 


(i) (triple solutions) when u{ui) > 0, there exist 0, 0 G M such that 

— u{(j)) (and n(0) ) takes three distinct values for 0 G (0, 0); 

— m(0) (and n(0) ) can be represented uniquely for 0 G (— 00 , 0) U (0, 00 ); 

— m(0 ) (and n(0)j takes two distinct values at cj) = 0,0; 

(^^) (uniqueness of monotone solutions) when (t{ui) < 0 , m(0) and n(0) can be 
represented uniquely for 0 G M. Moreover, m'(0) < 0 and n'(0) > 0 for cf) G M; 

(Hi) (unique and monotone solutions with inflection points) when a{ui) = 0 , 
m(0) and v{4>) can be represented uniquely for 0 G M. Furthermore, there exists 
0 G M such that 

— u\(j)) < 0 and n'(0) > 0 for 0 G (— C) 0 , 0) U (0, 00 ); 

— m'( 0 ) = = 00 at (j) = (j). 
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We are able to explain ci C 2 in the following manner. According to Theorem I3.21 the 
sign of a{ui) determines the properties, such as the number of solutions, of fl3.8jl . Here 
Ml is the largest positive root of p{u) = 0 and the discriminant as shown in fl3.5jl relies 
on di and gtj {i,j = 1,2). Consequently, Ui is determined once dt and gij {i,j = 1,2) are 
given. Now Cj and 6^ {i,j = 1,2) can be suitably chosen such that anyone of the three 
cases as described in Theorem 13.21 can occur. 


Theorem 3.3 (Existence of solutions to Poission equations with discontinu¬ 
ous nonlinearity). Assume that H C M"" is an open and bounded domain with smooth 
boundary dfl and 

(Al) p'{u) is piecewise continuous and discontinuous at finite points for u > 0; p'(s) ds] < 

oo for all C M; 

(A2) there exist constants K,L> 0, such that \p'{u)\ < L for \u\ < K and \p{u)\ > 
for |m| > K. 

Then the following Neumann problem 


Au = p'{u) in H, 

= 0 on dQ, 

has a unique weak solution u G 
Proof. For uq G M, let 

= \ / |VMpdx+ / {p{u) - p{uo)) dx. 

^ Jn Jo, 

As seen previously, it is sufficient to consider the minimizing problem 

min '0 (m). 
u£m{n) 

Using (A2), we have 

fi’iu) > 

> 

> 

> 


{p{u) - p{uo)) dx 


'\u\<K 


p{u) dx + 


'\u\>K 


-u^ dx — p{uo) |H| 


/ p{u) dx - p{uo) \n\ 

'\u\<K 

[ ([ p'(s) ds + p{0)\ dx - p{uo) \n\ 

'\u\<K \Jo J 

- (AL + |p(0)| + p(uo)) |H|. 


Therefore, there exists M G M such that 

M := inf fjiu) > —oo, 

ueH^{n) 


(3.9) 


(3.10) 


(3.11) 


(3.12) 

(3.13) 
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and we can find a sequence of functions {un{x)}^J^ G with lim^^oo'(/’(wn) = M. 

Now we show that ^jiu) is coercive, that is 

lini 4>{u) = oo, (3-14) 

|| li ||-|>00 

where ||m|| = ||M||L 2 (r 2 ) + ||VM||L 2 (o). Indeed, when ||VMn||L 2 (o) ^ C )0 as n —)■ oo, we clearly 
have fl3.14p . On the other hand, when ||VMn||L 2 (Q) is bounded and ||Mn||L 2 (Q) ^ oo as 
n —)■ oo, it follows that fl3.14p also holds. To see this. 


il{Un) > / {p{Un) - P(Mo)) dx 


/ p{Un) dx+ p{Un) dx - p{Uq) | 0 | 

'\Ur,\<K J\Ur,\>K 

1 r 

uldx - p(tio) |fi|. 


> -KL\a\ + - 


\Un\>K 


(3.15) 


Since | |Mn| |l2(q) —)■ oo as n —)■ oo and dx < |f2|, we conclude that dx —)■ 

oo and oo as n —)■ oo. The coercivity condition 03.1411 leads to the fact that 

there exists a constant Mq > 0 such that ||Mn|| < Mq (otherwise, the unboundness of 

||m„|| together with 03.1311 contradicts 03.14p i. It follows that the sequence {un{x)}'^^ 
is bounded in and there exists a subsequence {un-} of with 

Ujij u* in as j —?> oo. (3.16) 

We select a subsequence {unj^} of {un^} which converges in T^(fl), i.e. 

—!■ u* in L^(fl) as k ^ oo. (3-17) 

Since convergence implies pointwise convergence of a subsequence almost everywhere, 
we have 

Urij^ —^ u* a-e. in 11 as / —)■ oo, (3.18) 

where {unj^ } is a subsequence of {unj^}. We denote this subsequence obtained by 

To hnd a minimizer of our minimizing problem, it suffices to establish weak 
lower semicontinuity, i.e. 

i’iu*) < liminf'0(Mfc), (3.19) 

k—^oo 

which gives 


M = lim 'ijj{uk) 

k—^oo 

> lim inf 

k^oo 


= lim inf 

k^oo \ 2 


\'\/uk\‘^dx+ / {p{uk) - piuo)) dx 


-o \'^u*\^dx+ / {p{u*) - p{uo)) dx 
^ Jn Jn 

= ip{u*) > inf ipiu) = M. 
u€m{n) 


(3.20) 
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This shows that u* solves the minimizing problem fl3.11j) . Now we prove fl3.19p . Indeed, 

pU 

p{u) = / p'{s) ds + p(0) 

Jo 

is continuous since p'{s) G L^(r2*) for all hi* C M. Since Uk{x) —)■ u*{x) a.e. in hi as 
A; —oo, we have p{uk) —)■ p{u*) a.e.. Applying Patou’s lemma to obtain f^p(u*)dx = 
lim inffc^oo p(uk) dx < lim inffc^oo Jq p(uk) dx yields 


lim inf 'ip{uk) 

k^oo 


= lim inf I - 

fc—>oo V 2 




\S/uk\^dx+ / {p{uk) - p{uo)) dx 


In 


1 


>2 / j [p{u*) — p{uq)) dx = ip{u*), 


in 


in 


(3.21) 


where we have used the fact that | Vufcp dx —>■ | dx as k ^ oo. In fact, elliptic 
regularity of weak solutions ensure smoothness of u*. This completes the proof of the 
theorem. □ 


4 Trend to equilibrium 


In this section, the following initial conditions 


M(a;,0) = Uo{x) > 0, u(a:,0) = Vq{x) >0, a: G hi, (4.1) 

are imposed to determine completely the evolution. Since the hrst two equations of fll.3p 
are in divergence form, total charges are conserved: 


u= u{x,t)dx= / Uo{x)dx, v= v{x,t)dx= / Vo{x)dx. 


(4.2) 


The means Wi and W 2 of the charges u and v over hi are dehned by 


u ^ V 

= M- = M- 

respectively. The following two Lemmas are crucial in proving Theorem 14.31 We note 
that Csiszar-Kullback-Pinsker inequality (references refer to the ones cited in [3]) comes 
from information theory. 


Lemma 4.1 (Csiszar-Kullback-Pinsker inequality). For u,v G L^{Q), u,v>0, 

J log j — M + dx > Cx(fl) ^ y \u — v\dx'^ , (4.4) 

where Ck{^) is a constant. 
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Lemma 4.2 (Logarithmic Sobolev inequality ([!])). For y/u G u > 0, 


/ u log ^ dx + C{n, it) < Cl{^) / \'V\/u\'^dx, (4.5) 

Jn Jn 

where w = |^, u = J^udx, Cl{^) is the Logarithmic Sobolev constant, and C{n,7i) is a 
constant depending on n and tt. 


Consider the initial valne problem for fll.Sp with initial conditions u{x,0) = Uq{x), 
n(x, 0) = ^ 0 ( 3 ;)) 0(^)0) = 0o(^) and the Nenmann bonndary condition 


du dv dd) ^ 

^ LT ^ LT ^ X e t > 0. 
on ov ov 


(4.6) 


Under certain hypotheses on the initial conditions and the parameters appearing in 
jrsii. we show that the solution of the initial-boundary value problem fll.dp tends to wi 
and Wi in the sense. More precisely, we have 


Theorem 4.3 (Trend to equilibrium). Assume 


(■HI) %^(2 9 n - 9 i 2 - 921) > max ( - + 7i (2t9i + ^^ 2 )), O); 


(H2) ^(2g22 - 9 i 2 - 921 ) > max ( - |( 7 it 92 + 72 {‘^'^2 + ^9i)),0); 

(1-13) 0 < Hq < 00 , where Hq= (uq log ^ Uq log ) dx. 

Jn ^ U!l U!2'' 

In (HI) and (712), Cp is the Poincare constant. If a global-in-time solution {u{x, t),v{x, t)) 
of the time-varying problem (O, (SU), (USD exists, then it tends to the corresponding 
constant steady-state solution {wi,W 2 ) in the L^-norm as t ^ 00 , i.e. 

||m - WiIIlph), ||w - W 2 ||li(o)—^ 0, as t —)> 00 . (4.7) 

Furthermore, the time-varying solution {u{x,t),v{x,t)) converges in the L^-norm expo¬ 
nentially fast to its mean with explicit rate B: 

\\u - + \\v - W 2 \\l^n) <Ck~^ Hoe-^\ (4.8) 

where B will be specified in the proof, and Ck is the constant in the Csiszdr-Kullback- 


Pinsker ineguality (Lemma f.l). 

Proof. As in [5], we define the relative entropy functional H[u,v]{t) by 


f / U V \ 

Hlu, v]{t)= [u log ——\- V log — ) dx. 

Jn ^ wi W2I 


(4.9) 


By virtue of the elementary fact logz + ^ — 1 > 0 for z > 0, it can be shown that 
H[u,v]{t) > 0 for all f > 0. Indeed, 

0< /4log^ + ^-l)+nflog^ + ^-l)da; 

Wi U J \ W2 V / 


'n 


=H[u,v]{t )/ Widx 

Jn 

=H[u,v]{t). 


/ udx-\- W2dx — V dx 

In Jn Jn 


(4.10) 
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Note in particular that (u log dx, 4 (v log ^)dx >0 for all t > 0 . For simplicity 
of notation, let Xi = di Vu + m V 0 + gn u Vu + 5^12 u Vn and X2 = ^2 Vn + 792 V 0 + 
5^21 V Vu + g22 V Vn. We calculate the time derivative of H[u, n] to obtain 


dt H[u, v]= ut( log ^ + 1) +vt( log — + l)dx 
\ Wi J \ W 2 J 

= [ V-Xi flog—+ 1)+V-Xi flog—+ 1) dx 
Jn \ wi J \ W2 y 

= [ -Xi ■ vf log — + 1) -Xi • vf log — + 1) dx 
Jn \ wi y \ W2 y 

= - [ fxi- —+X2-—)dx 

Jn ^ u V y 

= — j ^di u~^ I Vm|^ + 'di Vcj) ■ Vu + gii | Vm|^ + gi2 Vn • Vu 

d2V~^ |Vn|^ + '^2 V 0 • Vn + 5^22 |Vn|^ + g2i Vu ■ Vnj dx 
= — di / M“^|VMpdx —d2 / n“^|Vnpdx 

t/ £7 t/ £1 

— 7i'di / M^dx —72'd2 / dx — (71'^2 + 72 di) / uvdx 

Jn Jn Jn 

-911 / |Vtip 91 - 922 / IVtipda: - (912 + 921) / Vu-Vvdx, 


(4.11) 


where we have used Green’s hrst identity and —A 0 = 71 u + 72 n to get V 0 • Vm dx = 
7 i /o dx+72 uv dx as well as V 0 -Vn dx = 72 dx+71 Jq uv dx. On the other 

hand, it is readily seen that f^uvdx = f^ {u—wi){v—W2) dx+J^ wi W2 dx. In particular, 
we have dx = Jq{u — wi)"^ dx + wf dx and v'^ dx = — W2Y dx + dx. 

As a consequence. 


dt Jdlu, v] = — di / u iVufdx —d2 / v |Vnfdx 

Jn Jn 

-71^1 / (u - dx - 72'd2 / {v - W2f dx 
Jn Jn 

-(71^2 + 72^1) / (n - - 102) dx 


5(11 / |VMpdx-5(22 / iVnpdx - (5(12 + 5'2i) / Vu • Vn dx - Wi), 
do do do 

( 4 . 12 ) 


where tci) = 71 •di tCi dx + 72 ^^2 Jq wl dx + (71 ■d2 + 72 di) wi m;2 dx. 

Under the hypotheses ("HI) and ( 7 ^ 2 ), applying the Cauchy inequality a b < lia^+b^), 
the Poincare inequality f^(u — Wi)^ dx < Cp iVu^dx, where Cp is the Poincare 
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constant, and the logarithmic Sobolev inequality in Lemma [4.21 we hnd 


dt H[u, v] < — 4 di C 


-1 


u 


u log — dx — 4:d2Cj 
Wl 


-1 


V log — dx 
W2 


- 7 i^i / {u - Wl) dx - ^2'd2 / {V-W2) dx 


In 


--( 71^2 + 72 ^^ 1 ) / {u-wi)'^dx--{'yi'd 2 + ^ 2 'di) / {v - W 2 f dx 
^ Jn ^ Jn 

- 9 n / \Vu\'^dx-g22 / 

Jn Jn 

+ ^{912 + 921) I \Vu\'^ dx+ ^{912 +921) I \Vv\^ dx - A{wi,wi) 

^ Jn ^ Jn 

<- 4 :mm{di,d 2 )Cj;^ H[u,v]-I-{'J21J1+71(2 191 +192)) [ {u - wi)‘^ dx 

^ Jn 

- 77(71^2+ 72 (21^2+ ^?i)) / {v-W2)‘^dx-A{wi,Wi) 

^ Jn 

+ 7t(5'i2 + 5'2i - 25(11) / iVnl^da;+ -(512+ 5-21-2522) / \Vv\‘^dx 
<-Amm{di,d 2 )C^^H[u,v]-h 72 'di+'yi{‘J' 9 i+' 92 )) [ {u - lUi)^ dx 


-(71^2+ 72 (2i 92+ ^?i)) / {v-W2)‘^dx-A{wi,Wi) + 

Z Jn 


Cp' 


Cp^ 


{912 + 921 - 2 9u) / (u - wii) dx + ^-(512 + 5'2i - 2522) / {V-W2) dx 


< — 4 min((ii, ^ 2 ) H[u, v] — A{wi, + 1 ). 

In fact, A{wi,wi) can be rewritten as 

A(wi,wi) =7ii9i +1 |n| +72'd2+2 1^1 + ( 71 '*^2 + 72191 )+ 1+2 |hl| 
I^|f 7 ii 9 i(^) + ( 71 ^ 2 + 72 i 9 i) ^+ 72 i 92 y 

\ \Wr>J Wn I 


(4.13) 


= tCn 


(4.14) 


Due to the compatibility condition (or the electroneutrality condition) of the Neumann 

problem for 0, we show that 7l(ihi, + 1 ) = 0. To this end, the boundary condition ^ = 0 
on dVl yields 


0 = 


'an 


d(j) 

du 


dx= V • (V0) dx = 7 i / udx + 72 vdx, 


(4.15) 


or 7 i lUi = — 721112 , which results in 71 ( 1111 , 1111 ) = 0. This leads to the following initial 
value problem involving a differential inequality 


{dtH[u,v\{t) < -B{di,d 2 ,CL) H[u,v]{t), t > 0, 
\ id[ii,ii](0) = ido, 


(4.16) 
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where B{di,d 2 , Cl) = 4 min((ii, ^ 2 ) C'l Hq = H[u, w](0) determined by (14.91) is a positive 
constant depending on uo{x) and vo{x). Gronwall’s inequality yields 


H[u,v]{t)<e-‘^^Ho, 


(4.17) 


where B = B{di,d 2 ,CL)- Thanks to (773), 0 < Hq < 00 . Since from the Csiszar- 
Kullback-Pinsker inequality (Lemma 14.ip . the L^-norm oi u — wi and v — W 2 can be 
controlled by H[u,v]{t). The same decay rate as in fl4.17p is given by 





(4.18) 


This completes the proof of the theorem. 


□ 


References 

[1] A. Arnold, P. Markowich, G. Toscani, and A. Unterreiter, On convex 
Sobolev inequalities and the rate of convergence to equilibrium for Fokker-Planck type 
equations, Gomm. Partial Differential Equations, 26 (2001), pp. 43-100. 

[2] V. Barcilon, D.-P. Chen, and R. S. Eisenberg, Ion flow through narrow 
membrane channels. II, SIAM J. Appl. Math., 52 (1992), pp. 1405-1425. 

[3] P. Biler and J. Dolbeault, Long time behavior of solutions to nernst-planck and 
debye-huckel drift-diffusion systems, Annales Henri Poincare, 1 (2000), pp. 461-472. 

[4] D. Chen and R. Eisenberg, Charges, currents, and potentials in ionic channels 
of one conformation. Biophysical journal, 64 (1993), pp. 1405-1421. 

[5] C. Conga, E. Espejo, and K. Vilghes, Remarks on the blowup and global 
existence for a two species chemotactic Keller-Segel system in European Journal 
of Applied Mathematics, 22 (2011), pp. 553-580. 

[6] B. Eisenberg, Y. Hyon, and C. Liu, Energy variational analysis of ions in water 
and channels: Field theory for primitive models of complex ionic fluids. The Journal 
of Chemical Physics, 133 (2010), p. 104104. 

[7] R. Eisenberg, From structure to function in open ionic channels. Journal of Mem¬ 
brane Biology, 171 (1999), pp. 1-24. 

[8] S. Fan, A new extracting formula and a new distinguishing means on the one vari¬ 
able cubic equation. Natural Science Journal of Hainan Teacheres College, 2 (1989), 
pp. 91-98. 

[9] C. Gardner, W. Nonner, and R. Eisenberg, Electrodiffusion model simulation 
of ionic channels: Id simulations. Journal of Computational Electronics, 3 (2004), 
pp. 25-31. 


20 





[10] T.-L. Horng, T.-C. Lin, C. Liu, and B. Eisenberg, Pnp equations with steric 
effects: a model of ion flow through channels, The Journal of Physical Chemistry B, 
116 (2012), pp. 11422-11441. 

[11] D. Horstmann, From 1970 until present: the Keller-Segel model in chemotaxis and 
its consequences. /, Jahresber. Deutsch. Math.-Verein., 105 (2003), pp. 103-165. 

[12] -, From 1970 until present: the Keller-Segel model in chemotaxis and its conse¬ 

quences. II, Jahresber. Deutsch. Math.-Verein., 106 (2004), pp. 51-69. 

[13] Y. Hyon, B. Eisenberg, and C. Liu, A mathematical model for the hard sphere 
repulsion in ionic solutions, Commun. Math. Sci., 9 (2011), pp. 459-475. 

[14] E. E. Keller and L. A. Segel, Traveling bands of chemotactic bacteria: a theo¬ 
retical analysis. Journal of Theoretical Biology, 30 (1971), pp. 235-248. 

[15] T.-C. Lin and B. Eisenberg, A new approach to the Lennard-Jones potential and 
a new model: FNF-steric equations, Commun. Math. Sci., 12 (2014), pp. 149-173. 

[16] -, Multiple solutions of steady-state Foisson-Nernst-Flanck equations with steric 

effects. Nonlinearity, 28 (2015), pp. 2053-2080. 

[17] B. Perthame, Transport equations in biology. Frontiers in Mathematics, Birkhauser 
Verlag, Basel, 2007. 

[18] T. Suzuki, Free energy and self-interacting particles. Progress in Nonlinear Differ¬ 
ential Equations and their Applications, 62, Birkhauser Boston Inc., Boston, MA, 
2005. 


21 




